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The internal structure of a composite fermion is investigated for a two dimensional parabolic 
quantum dot containing three electrons. A Yukawa screened Coulomb interaction is assumed, 
which allows us to discuss the evolution of the electron-vortex correlations from the Coulomb in- 
teraction limit to the contact potential limit. The vortex structure approaches the Laughlin limit 
non-monotonically through the formation of intermediate composite fermions in which a flip of the 
spatial orientation of the vortices with respect to the position of the electrons is observed. Only 
when we limit ourselves to the lowest Landau level (LLL) approximation the flip appears through 
the formation of an intermediate giant vortex at specific values of the screening length. Beyond 
the LLL approximation antivortices appear in the internal structure of the intermediate composite 
fermions which prevent the nucleation of giant vortices. We also studied the system of five electrons 
and show that the mechanism of the flip of the vortex orientation found for three-electron system 
is reproduced for higher number of electrons. 
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I. INTRODUCTION 

Theoretical interpretation of the fractional quantum 
Hall effectji (FQHE) observed at high magnetic field 
in the spin-polarized two-dimensional electron gas, is 
based on the properties of the Laughlin 2 . wave func- 
tion. FQHE for electrons is explained-, in terms of 
the integer quantum Hall effect for composite fermions, 
i.e. quasi-particles consisting of electrons with additional 
even number of bound vortices (or magnetic field fluxes). 
The vortices appear as zeros of the many-electron wave 
function when its phase changes by 2ir on a path around 
this zero. The electron in a composite fermion feels a 
reduced effective magnetic field as the bound vortices 
partly cancel the usual Aharonov-Bohm phase on a closed 
loop around the electronic The original problem consid- 
ered by Laughlin, 2 i.e. the diagonalization of the few- 
electron eigenequation in the basis of single-electron wave 
functions obtained in the symmetric gauge, is formally 
very similar to an electron system confined in a parabolic 
quantum dot. Only very recently wider attention was 
paid to the vortices in the quantum Hall regime of con- 
fined systems£i£iL2i2iiS and to the composite fermion the- 
ory for quantum dots^i*iSii^ In particular, the vortex 
distribution for Coulomb interacting electrons confined 
in quantum dots was investigated 5,6,7 using the exact di- 
agonalization technique and the reduced wave function 
imaging. The structure of vortices as obtained from such 
exact calculations differs significantly from the one as- 
sumed in the Laughlin wave function or in the composite 
fermion approach. It was found*^ that the vortices are 
not localized on the electron as assumed in the Laughlin 
state but stay in the neighborhood of electrons to which 
they are bound. On the other hand, Laughlin functions 
are the exact non-degenerate ground state wave functions 



for the case of short range interactions. Analytical proof 
of their exactness and uniqueness was provided** for po- 
tentials developed in series of V 2j (5 2 (r). The energy gap 
allowing the Laughlin liquid to be incompressible was 
identified* 5 - as due to the short-range component of the 
Coulomb interaction. 

The purpose of the present work is to investigate how 
the vortex structure is modified when the inter-electron 
interaction is taken from the Coulomb limit to the con- 
tact potential limit. We show that the vortices approach 
the Laughlin liquid limit in a non-monotonic fashion. 
Within the lowest Landau level (LLL) approximation for 
filling factors v < 1/3 intermediate composite fermion 
states are found with two additional vortices localized on 
the electron. Beyond the LLL approximation the inter- 
nal structure of the intermediate composite fermion turns 
out to be very complex with possible appearances of an- 
tivortices which can even be localized at the position of 
the electron. Within the LLL we found that only in the 
contact potential limit more than two extra vortices are 
localized at the electron position. 

In the present paper we focus our attention on the 
lowest number of electrons, i.e. N — 3, for which a 
nontrivial 5 internal composite fermion structure can be 
observed in the reduced wave function. Next, we verify 
the conclusions reached for N — 3 studying the vortex 
structure of a five electron system. To study the depen- 
dence of the structure on the range of the inter-electron 
interaction we assume that the electrons interact through 
a Yukawa potential 
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which in the large and small screening length (a) limits 
yields the Coulomb and the contact potential, respec- 
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tively. A potential of the form QJ is obtained for an ex- 
ternal Coulomb defect linearly screened by a three dimen- 
sional electron gas. 16 In fact the screening of the electron- 
electron interaction in electrostatic quantum dots results 
from charges induced on the metallic electrodes and is 
of a more complex form^i The screening of the electron- 
electron interaction by the image charges cuts off the long 
tail of the Coulomb potential. The contact potential limit 
corresponds then to the case of a negligible distance of 
the quantum dot to the metal gate in comparison to the 
dots size. 

This paper is organized as follows: Section II presents 
the theory behind the results which are given in Section 
III. Subsection III (a) contains the results calculated in 
the LLL approximation and the influence of the higher 
LL is described in subsection III (b) , results for five elec- 
trons are given in subsection III (c). Summary and con- 
clusions are provided in Section IV. 

II. THEORY 

The effective mass Hamiltonian of our system is 

H = E 2m* + Kxt(n) + E 

(2) 

where V e xt(r) = ^m*uu 2 r 2 is the parabolic confinement 
potential, and A is the vector potential. We adopt the 
GaAs effective mass m* = 0.067m e and dielectric con- 
stant e = 12.4. All the calculations were performed 
for Tluj = 1 m cV for which the oscillator length equals 
Iq = y/h/m*u> — 33.7 nm. The Schrodinger equation is 
solved using the exact diagonalization (ED) technique 18 
with the three-electron Slater determinants constructed 
from the single-electron Fock-Darwin orbitals^ We in- 
vestigated the ground-state magnetic-field induced angu- 
lar momentum and spin transitions of the three-electron 
system as function of the screening length in the presence 
of a perpendicular magnetic field [(0, 0, B) = V x A]. For 
a^oowe exactly reproduce the results of Ref.[2fi] (our 
parameters correspond to the interaction constant A = 
lo/a-B — 3.44, with as the donor Bohr radius). For finite 
values of a no interesting results are obtained: decreas- 
ing the screening length has the trivial effect of decreasing 
the strength of the interaction (A), the ground-state spin- 
orbital symmetry sequence remains unchanged, only the 
critical magnetic fields for the transitions between sub- 
sequent angular momentum states are shifted to higher 
values. 

We consider only the spin-polarized states of the magic 
angular momentum sequence 1 ^ [total angular momen- 
tum LTi being a multiple of 3h], which become ground 
states at high magnetic fields, after the maximum den- 
sity droplet decays. The results presented below were 
obtained mostly within the LLL (more precisely in the 
lowest Fock-Darwin bancU^ of zero radial quantum num- 
ber and nonnegative angular momentum) to keep a direct 



correspondence to the Laughlin wave function. In the 
discussion of the vortices we do not apply any magnetic 
field to the system without loss of generality for the wave 
function, since for a harmonic confinement potential the 
magnetic field simply rescales the electron coordinates of 
the wave function for a given LiAS 

*B#o(ri,r 2 ,r 3 ) = #5=0(71:1, 7r 2 , 7r 3 ), (3) 

with the scaling factor 7 = (1 + (lu c /2uj) 2 ) 1 ^ 4 , where 
lu c = eB/m* stands for the cyclotron frequency. Note, 
that property (3) implies that if, as generally accepted, 
the ground states at high magnetic fields are well approx- 
imated by the LLL, the approximation is not any worse at 
B = 0, where the high L states correspond to high excita- 
tions. Moreover, the eigenstates of the Hamiltonian writ- 
ten in the basis of Slater determinants built of LLL wave 
functions can be exactly identified with the eigenstates of 
the electron-electron interaction matrix operator. They 
are therefore the same for any constant A multiplying 
the interaction potential [Eq. (1)], even if for large A the 
LLL approximation can be arbitrarily bad^ In the cal- 
culations we consider screening lengths a > 0.1 nm. The 
delta-like interaction potential obtained for a — ► does 
not influence the energies or wave functions for a spin- 
polarized system because of the Pauli exclusion princi- 
ple. Consequently, for a — one obtains a multifold 
degenerate non-interacting ground-state. Since in the 
diagonalization these states (with very different vortex 
structure each) mix stochastically one cannot carry on 
the discussion of vortices for a screening length equal to 
zero. As a matter of fact, there is actually no need to take 
a strictly zero, since then the Laughlin function, as well 
as any other wave function constructed within the LLL, 
obviously corresponds to the degenerate ground state. 

A general form of the three-electron wave function in 
the LLL approximation is: 21 

z 2 , z 3 ) = ]T ViM 1 ^ exp (~ E ¥ ] 

j \ k=l / 

where A stands for the antisymmetrizer, Zk = Xk+iyk de- 
notes the complex, two-dimensional position of the fc-th 
particle, rjj are the linear variational parameters, j\, j'2, 
j'3 are nonnegative integers, of which not a pair is identi- 
cal, and ji +J2+J3 = L. The Laughlin wave function 2 for 
the angular momentum L = 3m (for odd m) is a product 
of the Jastrow factor and a Gaussian 

®i/ m {zi,Z2,z 3 ) = Y[(z k - z/) m exp (~ ^pf- ) . 

k<l \ n=l / 

(5) 

which is a special form of the general formula . In the 
Laughlin function the filling factor v = 1 /m is directly re- 
lated to the number of zeros m localized on each electron 
as well as to the angular momentum v = N(N — 1)/(2L) 
(for N — 3 electrons one has v = 3/L). Note that not 
all the states of the magic angular momentum sequence 
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FIG. 1: (color online) Positions of vortices of the conditional wave function calculated for the two electrons fixed at points 
(±Zo, 0) as function of the screening length a. Solid lines correspond to lower horizontal axis and show the positions of vortices 
on the x axis (y — 0). Dashed lines are plotted with respect to the upper horizontal axis and show the y coordinate of vortices 
localized on the x = —lo line. All the results were obtained in the lowest Landau level (LLL) with the exception of the blue 
curves plotted in (b) calculated beyond the LLL approximation with a fully convergent basis set. At the left side of (b) we 
show the electron-vortex orientation before (a > 0.44io) and after the formation of the giant vortex (a < 0.44Zo). The electron 
positions are marked with dots, and the vortices by crosses. 
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FIG. 2: (color online) Contour plots of the logarithm of the 
absolute value of the reduced wave function calculated for 
angular momentum L — 15 and the screening length a — Zo/2 
in the LLL approximation for the pinned electron positions 
(±Zo,0) (a), (±Zo/10,0) (b) and (=H + V 2 , 0) (c). Electron 
positions are marked by blue arrows in (a). 

can be represented by the Laughlin function, only those 
of odd L can. 

We investigate the zeros of the reduced wave 
function(£i£i£ constructed by fixing coordinates of two 
electrons z\ and Zi 

itzuzii 2 ) = *(zi,Z2,s), (6) 



where z is the test electron coordinate. The reduced 
Laughlin wave function is a complex polynomial of the 
test electron position (z) of degree 2m = |L, multiplied 
by a Gaussian. On the other hand, for a general LLL 
state (@J the reduced wave function is a complex poly- 
nomial of degree L — 1, resulting in more zeroes than 
occurring in the Laughlin wave function. The additional 
zeroes, commonly attributed to vortices bound to the 
test electron, are not localized close to the pinned elec- 
tron positions. Since one extra zero has to be attributed 
to the test electron itself, one obtains the total number of 
L vortices [for a general number of iV electrons the num- 
ber of vortices equals N/v = 2L/(N — 1)]. When higher 
Landau levels are included the reduced wave function 
depends also on the complex conjugate of the particle 
positions and larger exponent values in the polynomial 
appear, which increases the number of zeros and allows 
antivortices to appear m& 



III. RESULTS 

A. Lowest Landau level approximation 

Fig. 1 shows the position of zeros of the LLL reduced 
wave functions for two of the electrons pinned in the lo- 
cations (±Zrj,0) for states with L — 9, 12, 15, 18 and 21. 
In the presented range of x we focus only on the zeros 
located near the pinned electrons. In the Coulomb limit 
all the vortices are placed on the x axis£ For L = 9 [Fig. 
1(a)], as the screening length decreases, the two vortices 
bound to the electron approach its position, staying al- 
ways on the x axis. For a — 2 nm (a = O.O6Z0) the 
bound vortices are localized exactly at the electron posi- 
tion forming a giant vortex, characteristic for the Laugh- 
lin wave function. For L = 12 [see the black lines in Fig. 
1(b)] the giant vortex on the electron position is formed 
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earlier, i.e. for a — 0.44/q. However, for smaller screening 
lengths the two extra vortices leave the electron position 
(and the x axis) passing to the x = —l line [see the in- 
set in Fig. 1 (b)]. For still smaller a the vortices return 
to the electron position. A similar behavior is found for 
larger L. For states with L > 12, there are more than 2 
extra vortices bound to each electron and pairs of them 
collapse on the electron positions for specific L-dependent 
screening length values. Decreasing a beyond this value 
flips them to the x = ±Zo line approaching again the 
electron positions in the a = limit. Note, that the for- 
mation of the intermediate giant vortices is observed also 
for non-Laughlin states (even L) and that all these inter- 
mediate giant vortices have winding number three. For 
non-Laughlin states the number of zeros of the reduced 
wave function (L — 1) is odd, therefore a single vortex 
resides in the (0, 0) position in order not to break the 
symmetry. The position of this vortex for L — 12 and 
L = 18 is marked by the vertical line just to the left of 
the tick marks on the right hand side. We see that for 
even L, i.e. the non-Laughlin states, the number of vor- 
tices bound to the electrons is the same as in the closest 
Laughlin state with lower angular momentum (L — 3). 

The presented results are quite general in the sense 
that for N — 3 the vortex structure does not depend on 
the specific choice of the positions of the two fixed elec- 
trons in Eq. (6) but scales linearly as a function of the 
distance between the fixed electron positions, whether 
they are, or not, placed symmetrically with respect to 
the origin. This is demonstrated in Fig. 2, which shows 
a contour plot of the logarithm of the absolute value of 
the reduced wave function calculated in the LLL approx- 
imation for L = 15 and the screening length a = Zo/2. 
In Fig. 2(a) the two electrons are fixed at (±Zo,0) like 
in Fig. 1 (c). For a = Z /2 two of the bound vortices 
are localized perpendicular to the line between the elec- 
trons [cf. Fig. 1(c)] whose positions are marked with 
the blue arrows. In Fig. 2(b) the fixed electron coor- 
dinates were scaled down 10 times with respect to Fig. 
2(a), and in Fig. 2(c) the electrons were shifted to the 
left by Iq/2. The scalability of the vortex structure is 
evident from the form of the LLL wave function (4). If 
one decreases all the distances 7 times, one can take 7 L 
before the sum (since ji + j-2 + J3 = L) from the poly- 
nomial part, i.e. the one responsible for the appearance 
of the vortices. The invariance of the vortex structure 
with respect to the shift of the fixed electron positions is 
not evident from Eq. (4). In fact, the vortex structure 
of each of the Slater determinants is not invariant with 
respect to the shifts, and the invariance is only obtained 
in the entire basis containing all the LLL determinants. 
Since we are dealing with the harmonic oscillator poten- 
tial the exact wave function is separable into a product 
of the center of mass and relative motion wave functions 
* = F C M(z C m)G re i (zi ~z 2 ,zi-z 3 ,z 2 -z 3 ). The vortices 
are entirely due to the relative part. From the separable 
form it is clear that the vortices shift with the fixed elec- 
tron positions. It is quite remarkable that this feature 



of the exact solution is reproduced in the LLL approx- 
imation. For a general N the vortex structure remains 
invariant with respect to the size, position and orienta- 
tion of the polygon formed by the N — 1 fixed electrons 
as long as its shape is preserved. One cannot change the 
shape of the line segment linking the two fixed electrons 
for N = 3. But for N = 4 different vortex structures 
are obtained when the shape of the triangle formed by 
the fixed electrons is varied^ Fig. 2 shows also that only 
the vortex structure and not the reduced wave function 
scales with the positions of the fixed electrons. The non- 
scalability of the wave function results from the center of 
mass component of the wave function [or the Gaussian 
in Eq. (4)] 

The dashed lines in the upper part of Fig. 3 show the 
positions of the two remaining vortices, which did not fit 
into Fig. 1(a) for L — 9. These vortices are not bound 
to the electrons whose positions are pinned but belong to 
the test electron and disappear to infinity as the screening 
constant is decreased to zero. This behavior is expected 
since the number of vortices in the Coulomb problem is 
larger than in the limit of the Laughlin liquid (see the 
end of Section II). The red full curves close to the lower 
horizontal axis show the modulus of the corresponding 
reduced Laughlin wave function for y = with the elec- 
tron positions fixed at (±Zo,0). We see that the disap- 
pearing test-electron vortices are always localized beyond 
the region of the reduced Laughlin wave function local- 
ization. This is not always the case. Black solid lines in 
Fig. 3 show the position of vortices for the two electrons 
fixed at (— Zo, 0) and (— Zo/2,0). The outermost vortices 
are localized more closely to the electrons. We see that 
in this case, for decreasing screening lengths the vortices 
pass through the region in which the reduced Laughlin 
function (plotted in blue in Fig. 3) takes large values. 

In order to get an idea how well the electron-vortex 
correlations are described in the Laughlin wave function 
we project the reduced optimal wave function obtained 
within the LLL approximation to the one corresponding 
to the Laughlin many-particle wave function 

o _ <Vw i h(ju > m 

J zi,z 2 — = 1 (I J 

y/< >< # llZ2 IVi, Z2 > 

in which the positions of vortices as well as of the pinned 
electrons can be seen (tp L denotes the reduced Laugh- 
lin wave function). The overlaps calculated between ED 
wave functions of states with odd angular momentum and 
corresponding Laughlin wave functions are shown in Fig. 
0] In Fig. 0] (a) the two pinned electrons were placed at 
(±Zo, 0). Vortex positions for this case are shown in Figs. 
ID^a), 1(c), and 1(e). The overlap values increase mono- 
tonically for all the three Laughlin states with decreasing 
a, in the a — > limit they all achieve unity. Note also 
that the higher the angular momentum the smaller the 
overlap. In Fig. Q]we observe that for these three states 
the distance between the outermost bound vortex and 
the electron increases with L which is the reason why for 
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larger L the overlaps are smaller. Moreover, there are 
regions for L = 15 and L = 21 where vortices increase 
their distance from the electron with decreasing a, which 
is not reflected in the overlap plot, which apparently is 
more strongly determined by the decreasing distance of 
the electron from the outermost vortex. 

More interesting behavior is observed when the elec- 
trons are pinned closer to each other. We placed them 
in (— lo,0) and (— lo/2, 0) and the resulting overlaps are 
shown in Fig. 0] (b). For all three states there is a more 
sharp minimum as function of a. This is due to the ex- 
ternal vortex passing through the region in which the re- 
duced wave function is large as discussed in the context of 
Fig. [3J The minimal value of the overlap, which is almost 
zero, occurs exactly when the vortex position coincides 
with the Laughlin wave function maximum. The dis- 
tinctly different dependence of the overlaps for the fixed- 
electron positions is due to the fact that the reduced wave 
function is not scalable with the interelectron distances. 

The displacement of the vortex towards infinity for 
a — > and its effect on the reduced LLL wave function is 
illustrated in the contour plots of Fig. [3]for L = 15 with 
the fixed electron positions (— lo,0) and (— Zo/2,0). The 
position of the vortex is marked by a ★ in Figs. 5(a) and 
5(c). In Fig. 5(b), a corresponds to the overlap minimum 
[cf. Fig. 4(b)] the vortex is visible near x = 31q where 
it digs a hole in the wave function. Moreover, when the 
vortex is at the position of the Laughlin wave function 
maximum, it splits the LLL wave function into two al- 
most equal parts with opposite signs (see Fig. [fj). This 
makes the Laughlin function almost orthogonal to the 
ED wave function. When the vortex of the test electron 
passes beyond the maximum of the wave functions, the 
overlap starts to increase reaching unity for all the states, 
but this occurs earlier for smaller values of L. 

Another relevant quantity to be discussed is the pair 
correlation function (PCF), defined as 

W(z a , zt) =< #| 5 ^ - Z ") 6 ( Z J - z b)^ > ■ (8) 

PCF calculated for the Laughlin function JSJ) gives (up 
to a normalization constant) 

W L {z a ,z b ) = (z a - z b ) 2m exp[-(|z Q | 2 + \z b \ 2 )/l 2 Q ] 

x / dz(z a - zf m {z b - zf m exp(-H 2 /Z 2 ). 

Therefore, at small interelectron distances (z a — > z b ) the 
pair correlation function will asymptotically behave as 
W L (z a ,z b ) ~ (z a - z b ) 2m . We consider the PCFs with 
one particle fixed in \z a \ = 1.5Zo a s well as at the ori- 
gin z a =0. In the case of \z a \ — 1.5Zo we calculated the 
PCF for the other electron at the same distance from 
the origin (i.e. \z b \ — 1.5Z ) along an arc of 0.6/o length 
away from z a . Then, we fitted the results to a function 
of the form f(\z a — z b \) — a\z a — z b \ K . For the other 
considered pinned position [z a = 0) we repeated this 
procedure moving from the origin along a straight line 
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FIG. 3: (color online) Vortices and the conditional probability 
for L = 9 when the two electrons are pinned at (— lo,Q) and 
(—Zo/2,0) (black solid lines). The two outermost vortices for 
electrons in pinned at (±Zo, 0) are shown with dashed curves. 
The reduced Laughlin wave functions along the y = axis are 
also shown by red and blue lines, for the symmetrically and 
nonsymmetrically pinned electrons, respectively. 
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FIG. 4: Overlap integral of the conditional wave functions cal- 
culated for the lowest-energy state diagonalizing the Hamil- 
tonian in the LLL subspace and the state described by the 
Laughlin wave function JjJ. In (a) two of the electrons 
are fixed at positions (±Zo,0), and in (b) at (— lo,0) and 
Ho/2,0). 



of length 0.2Zq. The obtained results are shown in Figs. 
[7f a) and (b) . We found that the value of the exponent 
depends on the fixed electron position (z a ), which is not 
the case for the Laughlin wave function. For L = 9 [blue 
curves in Fig. Efa)] the fitted k value approaches 6 - the 
Laughlin limit - monotonically with decreasing a. For 
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FIG. 5: Contour plots of the absolute value of the reduced 
wave function for the state with L = 15 and the two elec- 
trons fixed at (— Zo,0) and (— io/2, 0) (indicated by crosses). 
(a)-(c) show the ED wave function as obtained within the 
LLL approximation for different values of a, (d) shows the 
Laughlin wave function corresponding to this state. One of 
the vortices bound to the test electron which crosses through 
the wave function's maximum is indicated by a star in (a) and 
(c), while in (b) it creates a distinct minimum at x = 31q. 




FIG. 6: Evolution of the absolute value of the reduced wave 
function at y — as function of a for the state with L = 15 
and two electrons pinned in (— lo,0) and (— Zo/2,0). 



L = 12, i.e. a non-Laughlin state, also the value of 6 is 
obtained in the a = limit. In this case three vortices 
become localized at the electron position [see Fig. 1(b)] 
like for L = 9. We also notice that the K fitted for differ- 
ent fixed electron positions [black solid and black dashed 
lines in Fig. 0a)] are both equal to 6 around 0.45Zo, 
i.e. when the intermediate giant vortex is formed [see 
Fig. 1(b)]. The intermediate giant vortex is therefore as- 
sociated with the appearance of a position independent 
k, which is characteristic of the Laughlin wave function. 
Note, that the k value calculated for z a = between the 
intermediate (a = 0A51q) and the final (a — 0) giant 
vortices possesses a local minimum. This minimum is 
related to those vortices which initially are moving away 
from the electron for a below the occurrence of the inter- 
mediate giant vortex (see the dashed lines in Fig. 1(b)]. 
However, for the k exponent calculated with \z a \ = 1.51q 
[black dashed line in Fig. 0a)] a maximum is observed 
below a = 0.45Zo- For L = 15 the intermediate giant 
vortex is formed around a = Iq [Fig. 1(c)]. The k values 
fitted for the two pinned electron positions approach one 
another near a = la [sec the solid and dashed red curves 
in Fig. 0a)], but the k value for \z a \ = 1.51q is larger 
than 6. As before, a more direct correspondence between 
the vortex positions and the fitted k value is obtained 
for the electron pinned at the origin. The loop that the 
two vortices perform in the (a, y) plane when the inter- 
mediate giant vortex decays into single vortices [see Fig. 
1(c)], has no effect on the k value for \z a \ = 1.5/o- On 
the other hand the loop is translated into a minimum of 
K as calculated for z a = [see the red solid line in Fig. 
0a)] . Both k values tend to the value of the Laughlin 
function, i.e. to 10, in the a — limit. This limit is 
achieved by k values calculated for L = 18 as well, since 
also here 5 vortices are found at the electron position 
in the contact potential limit [see Fig. 1(d)]. Again, for 
L = 18 the k value calculated for z a = is more sensitive 
to the actual vortex behavior. Local maximum (slightly 
above 6) is obtained [dashed curve in Fig. (b)] when 
the first intermediate giant vortex is formed [a ~ 1.61q, 
see Fig. 1(d)]. Another maximum is observed for the sec- 
ond intermediate giant vortex (a ~ 0.37/q)- The value is 
now considerably larger than 6, which can be explained 
by the presence of the other vortices localized in close 
proximity of the pinned electron. The third intermedi- 
ate giant vortex near O.l/o gives a plateau near k = 8.5, 
which then shoots up to 10, when the final giant vortex 
is formed. For L — 21 we observe again a local maximum 
in k calculated for z a = at a = 0.71q - an intermediate 
giant vortex position [see Fig. 1(e)]. For L — 21 we ac- 
tually do not observe the final giant vortex in the a — > 
limit, due to the problem of degeneracy of the ground- 
state as explained in Section II. However, the presented 
k values and the vortex positions for small a for which 
the ground-state is still nondegenerate, clearly indicate 
the giant vortex Laughlin asymptotic with seven vortices 
at the position of the electron. 



The close correspondence found between the PCF cal- 
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-1.206 


-0.813 


0.440 


- 


17 


61 


15.33831 


-1.246 


-0.779 


0.450 


0.0042 


19 


173 


15.33754 


-1.243 


-0.779 


0.421 


0.0073 


21 


392 


15.33732 


-1.238 


-0.783 


0.415 


0.011 


23 


761 


15.33722 


-1.231 


-0.786 


0.411 


0.014 


25 


1346 


15.33719 


-1.226 


-0.789 


0.409 


0.016 


27 


2213 


15.33717 


-1.222 


-0.792 


0.410 


0.017 


29 


3453 


15.33717 


-1.220 


-0.794 


0.413 


0.017 


31 


5158 


15.33716 


-1.223 


-0.795 


0.416 


0.016 
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TABLE I: Convergence of the results for L — 12 beyond the 
LLL approximation. E ni is the maximum energy of the non- 
interacting Slater determinants used for the construction of 
the basis set (B = 0). Second column lists the number of 
basis elements (K). E is the energy estimate for a = 1.48Zo- 
xi and x r are the positions of the bound vortices to the left 
and right of the electron fixed at the point (—lo, 0) as in Fig. 
1(b) for a = 1.48/o- a* is the screening length for which the 
distance between the vortices aligned in the horizontal and 
vertical directions is the same [see Fig. 8 (c)]. This distance 
(c) is listed in the last column. The first row of the Table are 
the results for the LLL approximation. Value for a* in the 
first row corresponds to the giant vortex. 



culated for z a — and the vortex behavior is quite re- 
markable. For L above the value for the maximum den- 
sity droplet, the charge density develops a minimum at 
the center of the quantum dot and the depth of the min- 
imum increases with L. Moreover, by fixing the posi- 
tion of one of the electrons at the origin one includes 
only those Slater determinants in which the zero angu- 
lar momentum Fock-Darwin state appears. The angular 
momenta of the two remaining orbitals must sum up to 
la + h = L (let l a < lb). From the asymptotic behav- 
ior of the single-electron orbitals at the origin [z l ] one 
should expect that the obtained re value is related to the 
lowest of all l a . Due to the applied fitting procedure one 
actually obtains a k > l a . For instance, the value of 6 ob- 
tained in the Laughlin limit for L = 9 indicates that the 
Slater determinants corresponding to angular momenta 
(0,1,8), and (0,2,7) do not contribute to the wave func- 
tion while (0,3,6) does. Expanding the Jastrow factor it 
is straightforward to check that the Laughlin wave func- 
tion contains admixtures of the (0,3,6) and (0,4,5) basis 
functions, but not of the (0,1,8) nor the (0,2,7) determi- 
nants. 



B. Beyond the lowest Landau level approximation 

In order to verify the calculated vortex structure in the 
neighborhood of the fixed electron we have performed ex- 
act calculations with a basis including higher Landau lev- 
els for L = 12. The basis was constructed in the following 
way. From all the Slater determinants built of the non- 
interacting Fock-Darwin states we picked up only those 
for which the energy at B — (see the discussion of the 
wave function scalability with the magnetic field given in 
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FIG. 7: (color online) Coefficient re determining the asymp- 
totic dependence of the pair correlation function W(z a , Zb) — 
\z a — z \ K for z a — * zt (see the text) for different total angu- 
lar momentum states. Solid curves were calculated as fits to 
the actual PCF values calculated on an arc of length 0.6/o for 
\z a \ = 1.5Zo, and the dashed curves on a line segment of length 
0.2Zq with one of the ends fixed at the origin for z a = 0. 



Section II) does not exceed a fixed energy value E ni . The 
number of basis elements K as function of E n i is listed in 
Table I together with the energy estimates obtained for 
an interacting system at a = 1.48Zrj. The first row of the 
Table corresponds to the LLL approximation. We obtain 
convergence of the energy estimate up to six significant 
digits. Fourth and fifth columns of the Table give the 
position of the vortices attached to the electron localized 
at point (— lo, 0) for the second electron pinned at (Zrj, 0) 
like in Fig. 1. The convergence of the position of the 
vortices is slower than the energy. Beyond the LLL ap- 
proximation for a = 1.48Zo and for a up to the Coulomb 
limit the distances between the electrons and the vor- 
tices are slightly larger than in the LLL approximation. 
The positions of vortices obtained with the most precise 
calculations are shown by the blue curves in Fig. 1(b). 



8 




FIG. 8: (color online) Contour plots of the logarithm of the absolute value of the L = 12 reduced wave function when the 
two electrons are pinned at (±Zo,0) calculated with a basis of 5158 Slater determinants including higher Landau levels. Plots 
(a-h) correspond to the screening lengths a/lo = 0.421,0.417,0.416,0.415, and 0.412, respectively. Red arrows point to the 
antivortex positions. 
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FIG. 9: (color online) Zoom of Fig. 1(b). Solid curves show 
the x-position coordinates of the vortices localized at the y = 
axis and the dashed lines are the {/-position coordinates of 
vortices localized at x = ±Zo line. Black curves correspond to 
the LLL approximation and the exact results are plotted in 
blue. 



Beyond the LLL approximation the wave function is 
nonanalytical and the exact number of nodes in the whole 
complex plane is not known a priori. However, we have 
found that within the range plotted in Fig. 1(b) the ex- 
tra nodes in the exact calculations appear only within the 
region where the LLL predicts the formation of the inter- 
mediate giant vortex. Fig. [8] shows the contour plots of 
the logarithm of the absolute value of the reduced wave 
function when the two electrons are pinned at (±Zo,0), 
for the range of a when the bound vortices flip their po- 
sitions from the x axis to the x = ±Zrj lines. Instead of 
the formation of the intermediate giant vortex, a state 
consisting of single separate vortices is formed. When 
the vortices approach the electron along the x axis [Figs. 
ISfa)], the node of the wave function associated with the 
electron elongates in the perpendicular direction and fi- 
nally splits into an antivortex localized at the electron 



position and two vortices localized at the x = ±Iq lines 
[Figs. IBfb)] placed symmetrically with respect to the 
electron position. For a certain screening length a = a* 
the distances between the vortices localized at the x axis 
and those localized at x = ±/ lines are equal [Fig. E£c)]. 
With decreasing a the vortices localized at the x-axis 
approach [Fig. IHId)] the pinned electron and annihilate 
with the antivortex localized therein. Eventually, we are 
left with a single vortex at the electron position and two 
vortices localized on the x = ±Iq line [Fig. |5^e)]. like in 
the LLL for a values such that we are between the in- 
termediate and the final giant vortices. This mechanism 
of the flip of the orientation of the vortices is found for 
all the wave functions calculated beyond the LLL with 
basis adopted according to the strategy explained above. 
Values of a* are listed in Table I. The corresponding dis- 
tances between the pairs of vortices (c) are given in the 
last column of the Table. Distance c initially increases 
with the size of the variational basis and finally saturates 
near 0.016Z . 

Fig. 9 presents a zoom of Fig. 1(b) for the range 
of a corresponding to the intermediate and final giant 
vortices. The blue curves are for the exact calculations. 
After the flip of the vortex orientation the results of the 
LLL and the exact calculations are nearly equal. Con- 
tribution of the higher LL becomes negligible when the 
electron-electron interaction is switched off. 



C. Five electrons 

The mechanism presented above for the flip of the vor- 
tex orientation is reproduced for higher number of elec- 
trons. To illustrate this we focused on the five-electron 
system at L = 35, i.e., a non-Laughlin state correspond- 
ing to a ground-state of the magic angular momentum 
sequence below the filling factor v < 1/3. This state is 
the counterpart of the L — 12 state for three electrons 
discussed in the context of Fig. 1(b). Calculations were 
performed in the LLL approximation. The plots of the 
logarithm of the absolute value of the reduced wave func- 
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tion are given in Fig. 10 for four electrons fixed at the 
corners of a square (±Z , ±^o)- Fig- 10(a) shows the case 
of the Coulomb potential, and Figs. 10(b-d) the case 
of the screened Coulomb interaction for a = 0.0889Zo, 
0.0643Z and 0.0222/ . In Figs. 10(b-d) we present the 
vortices near the electron localized at (Iq,Io). The vor- 
tices attached to the fixed electrons approach them along 
the diagonals of the square and form a giant vortex for 
a = 0. 0643?o [see Fig. 10(c)]. For smaller values of a 
the line along which the attached vortices are aligned 
is rotated over 90° as compared to Fig. 10(b) and is 
now perpendicular to the corresponding diagonal of the 
square [see Fig. 10(d)]. 




-2 -1 1 2 0.97 0.38 0.33 1 1.01 1.02 




FIG. 10: (color online) Contour plots of the logarithm of the 
absolute value of the reduced wave function calculated for five 
electrons at angular momentum L = 35 in the LLL approxi- 
mation. Plot (a) corresponds to the Coulomb interaction po- 
tential and plots (b,c,d) to the screening lengths a — 0.0889Zo 
a — 0.0643/o and a — 0.0222Zo- Positions of four electrons 
are pinned at the corners of the square (Ho, Ho)- Electron 
positions are marked by blue arrows. 



IV. SUMMARY AND CONCLUSIONS 

We have investigated the dependence of the vortex 
structure of a three-electron quantum dot on the range of 
the inter-electron potential. The Yukawa interaction po- 
tential can be changed continuously from the Coulomb 
to the contact potential Laughlin limit. The evolution 



towards the Laughlin liquid appears through the forma- 
tion of intermediate three- fold giant vortices at which the 
vortices flip their orientation with respect to the electron 
to which they are bound. In our discussion we relied on 
the reduced wave function where 2 electrons arc pinned 
and found that the screening lengths for which the gi- 
ant vortices are formed do not depend on the choice of 
the positions of the pinned electrons. Hence, for N = 3 
the giant vortices can only be created by manipulating 
the screening length and not the positions of the fixed 
electron. For N > 3 electrons the exact vortex structure 
in the reduced wave function depends on the shape of 
the polygon formed by the N — 1 fixed electrons. But 
the binding of the vortices to the fixed electrons for large 
L are independent of the exact location of the electron. 
It is the angular position of the bound vortices which is 
altered when we move the other fixed electrons. Nev- 
ertheless for N > 3, we find that the evolution to the 
Laughlin limit is also non-monotonic and is accompanied 
with flips of the vortex orientation and the formation of 
the intermediate composite fermion states. 

We found that the LLL approximation predicts the 
vortex positions quite accurately in the whole range of 
the screening length except for a values where the vor- 
tices approach closely the fixed electrons. For a certain 
value of the screening length we observe a flip of the vor- 
tex orientation. In general we found that this flip can 
be realized in four different ways: symmetry breaking, 
discontinuously, through a giant vortex or by the for- 
mation of antivortices. In the LLL approximation giant 
vortices (similar to the ones assumed in the Laughlin 
state) are observed at the orientation flip, even though 
vortices are expected to exhibit a repulsive behavior at 
close distances^ In the LLL approximation an antivortex 
can not appear because the number of zeroes of the re- 
duced wave function is fixed. When higher Landau levels 
are included extra vortices and an antivortex appear and 
annihilate preventing the formation of the giant vortex. 

The presented study of the pair-correlation function 
shows that the precise positions of the vortices with re- 
spect to the electrons are important for the physics of 
electron-electron correlations. The number of bound vor- 
tices in the close neighborhood of the electron is trans- 
lated into an asymptotic power-law form for the pair 
correlation function around the pinned electron position. 
For the giant vortices, i.e. for the intermediate composite 
vortex states, the electron-electron correlations acquire 
properties similar to the ones described by the Laughlin 
wave function. 
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